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Abstract 

Chcrn insulators are periodic band insulators with the property that 
their projector onto the occupied bands have non-zero Chern number. 
Chcrn insulator with a honiogeneoTis boundary display continuum spec- 
trum that fills the entire insulating gap. The local density of states cor- 
responding to this part of the spectrum is localized near the boundary, 
hence the terminology edge spectriim. An interesting question arises, 
namely, if a rough boundary, which can be seen as a strong random 
potential acting on these quasi 1-dimensional states, would destroy the 
continuum edge spectrum. This paper shows how such question can be 
answered via a newly formulated abstract framework in which the ex- 
pectation value of the current of a general observable is connected to the 
index of a specific Fredholm operator. For the present application, we will 
connect the expectation value of the charge edge current with the index of 
a Fredholm operator that remains invariant under arbitrary deformations 
of the boundary. 



1 Introduction 

For magnetic Schrodinger operators in half-plane, Hatsugai established a 
fundamental result [1, 2] that says that the number of conducting channels 
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forming in the gap of the bulk system due to the presence of an edge is 
equal to total Chern number of the bands below that gap. The work by 
Hatsugai applies only for rational magnetic fluxes and homogeneous edges 
with Dirichlet boundary conditions. Almost 10 years later, using an advance 
mathematical machinery, Kellendonk, Richter and Schulz-Baldes established 
[3] a new link between the bulk and edge theory, which ultimately allowed 
them to generalize Hatsugai's statement to half-plane magnetic Schrodinger 
operators with weak random potentials, irrational magnetic fluxes and gen- 
eral boundary conditions (the result was first announced in Rcf. [4]). The 
result was later extended to continuous magnetic Schrodinger operators in 
Ref. [5]. These breakthrough results are especially important since they 
directly relate the quantization of the edge currents to a new topological 
invariant, the index of a specific Fredholm operator. Using this new invari- 
ant, one can directly explore the topology of the edge states under quite 
general (and physically relevant) boundary conditions, without the need of 
any gedanken experiments that employ artificial boundary conditions. 

Under similar assumptions, the equality between bulk and edge Hall con- 
ductance was also demonstrated by Elbau and Graf, soon after the publi- 
cation of Ref. [6], this time using more traditional methods. In Rcf. [7], 
which treats continuous magnetic Schrodinger operators, the edge appears 
at the separation between a left and a right potential, potentials that can as- 
sume quite general forms, in particular, they can include strong disorder. A 
similar result was established for discrete Schrodinger operators in Ref. [8]. 
We should mentioned that certain regularization of the edge conductance 
was necessary for the case of strong disorder. It is also important to notice 
that these papers [6, 7, 8] give a direct relation between the edge and bulk 
conductance and they do not discuss and edge index. However, in many 
situations, the edge index introduced in Ref. [5] is quite useful [13] and this 
is precisely what we will discuss in this paper. 

The works we mentioned above studied the Integer Hall effect for simple 
magnetic Shcrodinger operators in half-plane. In 1988, Haldane introduced 
a new model that exhibits the Integer Quantum Hall effect even thought 
the net magnetic flux per primitive cell is zero [9] . This strange property is 
typically present in tight-binding models with more than one quantum state 
per site and with broken time-reversal symmetry. It is a consequence of the 
non-trivial Chern number associated to the occupied bands of the model. 
This work lead to the discovery of a new class of insulators that now bear 
the name Chern insulators. It was argued in the literature [10, 11] that, due 
to their non-trivial topological properties, finite samples of Chern insulators 
must posses edge conducting channels that wind around the sample similar 
to what happens in the magnetic case. Thus, Chern insulators will display 
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remarkable edge physics, without the need of a macroscopic magnetic field. 
We should mention that no crystals of Chern insulators were discovered so 
far. It was argued, however, that by patterning ferroelectric materials one 
could fabricate a photonic crystal [12] with similar properties. 

Theoretically, Chern insulators are as ubiquitous as the normal insulator. 
Generally, thight-binding Hamiltonians with more than one quantum state 
per site and with broken time reversal symmetry can be always tuned to 
exhibit bands with non-zero Chern number. One of the simplest example of 
a Chern insulator is given by 

H= E {|n,l)(n + b,,2| + |n + bi,2)(n,l|} 

n,i=0,l,2 

(1) + E {t|n,l)(n + bi,l| + t*|n + bi,l)(n,l| 

n,i=l,2 

+t*\n, 2) (n + hi, 2| + t\n + hi, 2) (n, 2|}, 

where n labels the sites of a planar lattice generated by bi and b2, bo=0 and 
Im[i]7^0. This model is unitarily equivalent to the original Haldane's model 
and, in the bulk, it exhibits two bands with Chern numbers equal to ±1. 
The edge structure of the model was studied in great detail in Ref. [11] and 
the numerical calculations show one edge channel forming along the edge 
of the sample, in accordance to what one will predict from the topology of 
the bulk bands (in these calculations one rather sees two channels since the 
calculations were performed for slabs rather than half-plane samples). The 
calculations were, of course, performed for homogeneous edges. Our main 
question here is what happens to the edge channels for a rough boundary. 
The question is interesting since the edge channels are localized near the 
boundary and the rough boundary can be seen as a strong random potential 
acting on these quasi-onedimensional states. 

While the question is interesting in itself, by solving it we want to demon- 
strate an application of the abstract machinery given in Ref. [13], where we 
developed an axiomatic framework for current quantization. In the process, 
we identified the minimal conditions in which such quantization takes place. 
We claim that this minimal conditions bring significant simplifications over 
the previous works and the main goal of the present paper is to show how 
standard it is to demonstrate that these minimal conditions are met in Chern 
insulators. For example, in Ref. [5], the authors rely heavily on the Gaussian 
decay of the heat kernel for bulk magnetic Schrodinger operators and the 
estimates throughout the paper are quite involved. Here we will show that 
just simple exponential decay estimates on the bulk resolvent will do the job. 
Moreover, we will demonstrate that, once these exponential decay estimates 
are established, the minimal conditions stated in Ref. [13] can be derived by 
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Figure 1 . The non-zero hopping terms of Hq link site n with 
sites that are confined in the square drawn with the dashed 
line. The square contains (2L + 1) x (2L + 1) number of sites. 



following several standard, model independent steps. Thus, the present ap- 
plication completes the abstract formalism of Ref. [13] by showing that the 
required minimal conditions can be demonstrated via a somewhat universal 
technique. We hope that the present simplifications will make the technique 
more accessible to the theoretical physicists and that the generality of the 
proofs will encourages a search for new systems where the technique could 
lead to non-trivial insight. 



2 Definition of the physical systems 

In this section we define a bulk system, whose dynamics is determined by a 
very general tight-binding Hamiltonian Hq. Then we define and parametrize 
what we call a rough boundary T. A probability measure is introduced 
over the set of all admissible boundaries, which will be used for averaging. 
An auxiliary Hamiltonian H and the Hamiltonian H^ for the system with 
a boundary T is also introduced in this Section. We describe the most 
important properties of these systems, whose complete derivations will be 
given in the Appendix. 
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2.1 The bulk system. 



We consider a tight-binding Haniiltonian Hq defined on the 2-dimensional 
lattice Z^, which has K number of quantum states per site. We use n=(ni, 712) 
to represent the sites of the lattice. The Hilbert space H is spanned by the 
vectors 

(2) |n,a),nG Z^,a = 1, . . . , K; {n,a\n' ,a') = Sn,n'Saa'- 
We adopt a very general form of the tight-binding Hamiltonian: 

K 

(3) ^0= 1; E E [C^|n,a)(n + m,/3|+C;|n + m,/?)(n,a|]. 

neZ2 m a,/3=l 

The sum over m runs over a finite number of points of Z^. To be specific, 
we consider that these points are contained in a square of (2L + 1) x {2L + 1) 
number of sites, centered at the origin (sec Fig. 1). An important assumption 
on Hq is the existence of an insulating gap A=[E-,E^] in the spectrum. 
Under these conditions, the bulk system has the following general property. 



Proposition 1. The resolvent of Hq has the exponential decay property: 

(4) {n,a\{Ho-zr^\n',l3) < ^ " ' 

dist(z,cr(ifo)) - C[q) 

where (^(q) is specified in Eq. 65. z is an arbitrary point in the complex 
energy plane, away from the spectrum a{Ho) of Hq. The above estimate 
holds only when q is small enough so that the the denominator in the right 
hand side is positive. 



We defer the technical proof to the Appendix. But we want to mention 
here that the exponential decay or localization of the resolvent away from 
the spectrum of the Hamiltonian is a very general property, which can be 
demonstrated through a variety of the techniques. For our concrete exam- 
ple, the property follows from a simple application of the Combes-Thomas 
principle [15]. The same principle can be used to obtain sharp exponen- 
tial decay estimates for continuous Schrodinger operators in the presence of 
magnetic fields [16]. 



2.2 The system with a rough boundary 

Let us consider a line F like in Fig. 2, where the only constraint we impose is 
that all points of F be confined within a distance D from the line ni=0. The 
line F can be described by a sequence {jnjm where 7„ gives the deviation 
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Figure 2. Exemplification of a rough edge T and of the def- 
inition of 7„. 



of r from the axis ni=0 at the row 712 =n of our lattice, as exemplified in 
Fig. 2. We have 7n £ I, with I = {-D + 1/2, -D + 3/2,..., D - 1/2}. 
We recall that we use n=(ni, to describe a point of the lattice. Thus, F 
can be considered as a point of the set $1 = F={. . . , 7_i, 70, 71, • • •}. 

On the set O, we introduce the product probability measure, denoted by 
dT, which is the infinite product of the simplest probability measure z/ on 
n) — X^ngx/C^)] /(^) being any function defined on I. We 
remark that dT obtained this way is ergodic relative to the discrete transla- 
tions along the vertical direction of our lattice. We will use the probability 
measure dT to average over all possible contours F. 

We assign a sign to the points of the lattice, depending on their position 
relative to the curve F: 

, . _ J —1 if n is to the left of F 

^ ' - I +1 if n is to the right of F 

The Hilbert space Ti decomposes in a direct sum 7i = 7i+ © Ti- , where 

(6) Ti.± = span{\n, a),Sn = ±1} 

We first introduce a very useful Hamiltonian H, which is obtained from 
Hq by erasing the hopping terms that cross the contour F. Using the sign 
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function introduced above, we can write the Hamiltonian H as: 

(7) H= Yl Y.^r^p\'^,a){n + m,P\ + T^^\n + m,P){n,a\]. 

Sn Sn+m>0 a/3 

It is also important to notice that H = Hq — AF with: 

(8) /W= E[ra^l"'")(" + "^'/^l+ra;|n + m,/3)(n,a|]. 

<0 al3 

The sums involve only sites that are localized near the edge, more precisely 
those with the first coordinate ni in the interval [-Z)-L+1,D+L-1]. It is also 
interesting to notice that H describes two decoupled systems with edges. 
Indeed, H decomposes into a direct sum H=H-®H^, where 

H- : T~L- — > Ti-i 

(9) i?_= E E[rS>|n,a)(n + m,/3|+r™;|n + m,/3)(n,«|] 

<0 ajS 

and 

H+ ■ 'H+ —>■ H+, 

(10) H+= E E[r™^|n,«)(n + m,/?|+r™;|n + m,/3)(n,a|]. 

>0 a/3 

The system with the boundary T is defined by the Hamiltonian Hy=H^. 
We list three important properties of the system with the boundary. 

Proposition 2. Let (j){e) be a smooth function with support in the spectral 
gap of Hq. We consider that the support of (j) is separated from the spectrum 
of Hq by at least a distance 6>0; 6 will be considered fixed from now on. 

1) Consider two lattice points, n and n' in the "+" side of the lattice. 
Then there exists A{q)>0, independent of n and n' and T, such that: 

(11) \{n,a\cPiHr)\n',(3) < Aiq)e-'i^^^+<\ 

for q small enough. 

2) Fix a positive integer N. There exists Bi\r>0, independent of n and 
n' and F, such that: 

(12) |(n,«|</)(//r)|n',/3) < _ ^.[^ ^^at ■ 

3) There exists CN{q) > 0, independent of n and n' and F, such that: 

„-g(ni+n'i) 

(13) \{n,a\4>{Hr)\n',P)\<CNiq)- 



|n2-n'2| + l)^' 
where q must be taken small enough. 
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We defer the technical proof to the Appendix. Here we want to make the 
following comments. The first property states that 0(i^r) is exponentially 
localized near the boundary. This property follows from the simple identity: 



(14) R{z) = Ro{z) + Ro{z)AVRo{z) + Ro{z) AV R{z) AV Ro{z) , 



where Rq{z)={Ho—z)~^ and R{z)={H — z)~^ . The identity written in Eq. 14 
reveals the following simple principle: if the resolvent of the bulk Hamilton- 
ian has the exponential decay property, then the difference R{z) — Ro{z) is 
exponentially localized near the edge, which is assured by the fact that AV 
is strictly localized near the boundary. Then the exponential localization of 
4>{H-p) follows from Stone's formula: 



where the inclusion of Rq has no effect, since the integral over A is con- 
strained inside the spectral gap of Hq. We point out here that i?(A + ie) 
is evaluated infinitely closed to the spectrum which can be singular. The 
Appendix shows how to deal with this problem in a very general way. 

The second property follows from an estimate on H similar to that of 
Proposition 1 and from the rapid decay to zero of ({q) as q goes to zero. 
The third property is a direct consequence of the previous two properties. 
As we shall see, these three properties assure that the minimal conditions 
stated in Ref. [13] are satisfied by the Chern insulators. 

3 The problem and the result. 

We define first the central observable, the operator yr defined on Hr, which 
gives the vertical coordinate: 

(16) yr\n,a) = n2\n,a), n = {ni,n2). 

The index T is there to indicate that the operator is defined on Ti-p. The 
observable yr has discrete spectrum, a{yr)=Z, each eigenvalue being infin- 
itely degenerate. Let TTr{n) denote the spectral projector of yr onto the 
eigenvalue n. We have the following explicit expression: 



(15) (f){H) = lim / Im[i?(A + ie) - Rq{X + ie)](^(A) — , 




R 



(17) 




ni>7n,« 



where {'jn}nez is the sequence corresponding to the contour T, as discussed 
in the previous Section. 
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On the large Hilbert space H, we can implement the discrete lattice trans- 
lations group along the vertical direction by: 

(18) Un\{ni,n2),a) = |(ni,n2 - n),a). 

The discrete lattice translations along the vertical direction can also be 
extended to a map t„ acting on the space 0, of all possible contours T. The 
map tn simply shifts a contour downwards by n sites. 

Let us collect now the important facts into the following list: 



• We have defined the family of self-adjoint Hamiltonians Hr'.Hr^'Hr, 

Ten. 

• The observable y obeys: 

(19) UnyrK = ytnV + n 

• The family of Hamiltonians Hy is covariant, namely, it„ is an isometry 
that sends Hy into HtnT and UnHj-u'^ = -f/t„r- 

• On the set 0, we defined the probability measure dT, which is ergodic 
and invariant relative to the mappings tn- 



We now define the trace (notation tro) over the states of zero expectation 
value for yr'. 

(20) tro{A} = TV{7rr(0)A7rr(0)}, 
and we use tro to define the current: 

(21) Jr = tTo!^p{Hr)^^^=itTo{p{Hr)[Hr,yr]}. 

Here p{e) is the statistical distribution of the quantum states. Since we are 
interested in the contribution to the current coming from the edge states, 
we assume that the support of p(e) is entirely contained in the interval 
[E^ + S,E^ — 5] and that J p{e)d€=l. Wc assume that p(e) is a smooth 
function. Note that tro above is finite precisely because of the properties 
stated in Proposition 2. 



Main Theorem. Let F{e) = p{e). Note that F(e) is smooth and equal 
to 1 below E- + S and to above E^ — 6; also F' (e)=— p{e) . Wc recall that 
E± are the edges of the spectral gap of the bulk Hamiltionian Hq. Using 
the spectral calculus, we define the following unitary operators: 

(22) Ur = e-^'''^^"^\ 

If TTjt is the projector onto the non-negative spectrum of yr, then: 
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Figure 3. Example of a finite deformation of contours: F <-> F'. 

(23) / dT Jt = ^ Index {7r+C/rvr+} . 

The Index is an integer number, independent of the shape of /9(e) or the 
contour T. 



The proof of the statement is given in the following sections. We continue 
here with a discussion of the Index, more precisely its invariance property 
and how to compute it. The Index is defined on the class of Fredholm 
operators as: 

(24) Ind^ = dim Ker[A] - dim Ker[A*], 

for A : X ^ Y, with X and Y two Hilbert spaces. The Index has several 
important properties: 

1) Ind^* = -IndA. 

2) IndBA = IndB + IndA, ior A : X ^ Y and B : Y ^ Z two Fredholm 
operators. 

3) Ind(^ + C) = IndA if C is compact (in particular, if C is finite rank). 

4) The Index is invariant to norm-continuous deformations of the opera- 
tors inside the Fredholm class. 



Based on these general properties, we argue that the Index is independent 
of the contour F. We consider a finite deformation as shown in Fig. 3, where 
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the contour T is deformed to the left, in a finite region, into the contour T' . 
We can equivalcntly say that T' is deformed to the right into T. In any case, 
Hr C Hr'- Let i : Hr Ht' be the inclusion map and p : Hr' — > Ti-r be 
the projection 

(25) p an\n,a){na\= an\n, a){na\. 

Note that p = i*, po i = 1-^^, and io p = tt-^^, (the projector from Hr' to 
TCr)- Because the two curves differ only in a finite region, dim i<'er[p] < oo 
(i^er[i] is empty) and the two operators are Fredholm. The insertion of 
TTp [/pTTp in the Hr' given by z o vTp t/rVTp o p has same Index since 

(26) Ind i o -K^Urir^ o p = Ind i + Ind vrjt C/rvrr + Ind i* 

and the first Index and last Index cancel each other. The action of i o 
TTpt/rTTp op on the Hilbert space Hr' is very simple, 

i O TTp JTrTTp op = (i O TTp O p)(^i O [/p O p) (i O TTp O p) 
^^"^^ = (i O TTj!: O _p)e-2TiF(ioi?rop) (^ioTT+Op). 

We notice first that the projectors iott^ op and ttf' differ by a finite rank 
operator. Consequently, 

(28) Ind {ioTT+o p^e-^^iF{ioHrop) o o 
and 

(29) Ind Tr+e-^'''^^'°"^°PK+ 

are the same. We also notice that i o Hr o p, acting on Hr', has the same 
expression as Hr, thus i o Hr o p and H^' differ by a finite set of hopping 
terms. As the estimates of the next Sections will show, we can continuously 
switch on the missing hopping terms in Hr until it becomes identical to Hr' 
and keep 

(30) 7rit,e-2''^^(^deformed)7r+ 

in the Fredholm class during the process. Due to the invariance of the Index 
under norm-continuous deformations, we can conclude 

(31) Ind TTrUr'^r = Ind irr'Ur'TTr'- 

Using successive deformations of the type describe above, we can always 
deform a contour into another and the conclusion is that the index is inde- 
pendent of r. 

The Index can now be computed by taking the contour F as a straight 
vertical line, in which case we can use the translational invariance, more 
precisely, the Bloch decomposition. Let us denote by Fq such a straight 
contour. Definitely the theorem stated above applies equally well to the 
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Figure 4. The figure lists all possible ways an edge band 
(blue lines) can cross the energy window [E^ + 6, -E+ — 5] . 
Below each configuration the figure presents the value of the 
coefficient 



case when the set O reduces to one point, the contour Tq (all we have to do 
is to take D=0). Then we have the following practical way of computing 
the Index: 

Ind7r+ C/roVrro = Tr{^ro(0)p(i?ro)[^^ro, yro]7rro(0)} 

TT 

= Y1 I P{<^n,k)dken,k dk = YjOln. 
n k=-Tr n 

where e„^fc are the edge bands and q„ are defined in Fig. 4 (we used here 
the fact that / p{e)de=\). In other words, the Index gives the difference 
between the number of forward and backward moving edge bands. This is 
a very simple but important result because the edge bands can be easily 
computed numerically for the homogenous case. Moreover, if necessary, one 
can further simplify the computation of the Index by considering continuous 
deformations of the bulk Hamiltonian that keep the insulating gap opened. 
Using the fundamental result of Kellendonk, Richter and Schulz-Baldes [3] , 
we now know that the number computed in Eq. 32 is equal to the Chern 
number of the bulk states below the insulating gap. 

As a final remark for this Section, note that our main statement is about 
the average of the edge current and not the current itself. However, since 
the family {ffrjreo is covariant relative to translations, which act ergodi- 
cally on Q, the spectrum of {-ffr} is non-random. This implies that, if the 
edge spectrum becomes localized for a non-zero measure set of it will 
be localized for all contours, except a possible zero measure set. But this 
cannot happen, exactly because the average of the edge current is non-zero 
for Chern insulators. This allows us to conclude that the rough edge cannot 
destroy the edge conducting channels. 



(32) 
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4 Proof of the Main Statement 

Before starting the proof of our main theorem, we collect a set of technical 
facts, which all follow from the estimates given in Propositions 1 and 2. 

4.1 Technical results 

Along this paper, the following notations || || and || ||hs represent the operator 
and the Hilbert-Schmidt norms, respectively. 

Proposition 3. The following statements are true: 

1) {Ut — /)7rr(n) are Hilbert-Schmidt. Moreover, their Hilbert-Schmidt 
norm is less than an upper bound, which is independent of T. 

2) With the notation 

(33) Kr(n,n') = ||7rr(n)(C/r - /)7rr(n')||L, 
there exists Gjv > 0, independent of F, sucht that 

(34) KrM < j^^^^, 
for any positive integer N. 

3) If TTp denotes the projector onto the negative spectrum of yr and 
Sp = TTp — TTp , then [Sr!?7r] are Hilbert-Schmidt. Moreover, their 
Hilbert-Schmidt norm is less than an upper bound, which is indepen- 
dent of r. 

4) [yr, Ur]Trr{n) are Hilbert-Schmidt. Moreover, their Hilbert-Schmidt 
norm is less than an upper bound, which is independent of T. 

Let us also state a fundamental property of tro, essential for the proof of 
our main statement. 

Proposition 4. If {^rjren and {i?r}ref2 are two covariant families of 
operators such that 7rr(0)Ar, Ar7rr(0), 7rr(0)i?r and -Brvrr(O) are Hilbert- 
Schmidt with uniformly bounded norms, then: 



(35) J dr tro{ylrSr} = J dT tTo{BrAr} 



< oo. 
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Figure 5. Illustration of several definitions used in the text. 



We also need to introduce an approximate spectral projector onto the 
support of yo(e). For this, we consider a smooth function G(e) which is equal 
to below E- + (5/2, to 1/2 on the interval \E^ + 5, Ej^ — and to 1 above 
Ej^ — 6/2 (see Fig. 5). In this case, 

(36) = ^ - ^-'^'''^''^^) 

leaves invariant the states corresponding to the spectral interval of support 
of /o(e). By construction, tt^ has similar properties as Ut — I, which are listed 
stated below. Note that tTs depends on F, in fact {vTsjr form a covariant 
family. The dependence on F will not be stated explicitly. 



4.2 The Proof 



The proof makes use of the following non-commutative version of the Residue 
Theorem, whose proof can be found in Ref. [13]: 

A Non-Commutative Residue Theorem. Let f(z) be analytic in a strip 
around the unit circle. If {Ur}ren a covariant family of unitary operators 
such that {Ur — /)7rr(0) and [yr, f^r]7rr(0) are Hilbert-Schmidt, then: 

(37) JdT tro{(/(t/r) - /(/))[yr, f/r]} = h J dT tio{{U^ - I)[yr,Ur]}, 
where bi is the coefficient appearing in the Laurent expansion: 

(38) f{z)= ^ ^ a„z^ 

n=l,oo n=0,oo 



Now from Proposition 4 we know that [Sp, Ur] is Hilbert-Schmidt. Then 
TTp JTr^Tp is Fredholm and we can use Connes' result in non-commutative 
geometry: 

(39) Ind{7r+?7r7r+} = -^Tr{Sr[Sr, ?7f^] pr, f/r]}, 
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which, by using elementary manipulations [13], can be reformulated in a 
slightly different format: 

(40) Ind{7r+C/r7r+} = - ^ ^{4iUr - ^)Pr, Ur]4}. 
We use the projectors '7rr(n) to expand 

(41) Ind{7r+f/r7r+} = El^K(n)(i7^ - /)[Er, C/r]7rr(n)}. 

n 

We now consider the average over V. Since the index is independent of 
r, the average over T can be omitted for the left side. On the right hand 
side, we use the fact that the trace of trace-class operators is invariant to 
unitary transformations and the fact that the measure dT is invariant to the 
mappings t„, to write: 

IndlTrjtC/rTrf} 
= -i/drElV{w(n)(C/^ -/)[Sr,C/r]7rr(n)<} 

n 

(^2) = / dr E TV{7rt„r(0)(C/;„r " /)KSr<, C/i„r]vri„r(0)} 

n 

= -I /drElV{7rr(0)(C/f - /)KEt_„r<, Ur^m 

n 

The key observation at this step is that: 

(43) ^n„Si_„r< = 2yr + /, 

n 

whose graphical representation is given in Fig. 6, which leads to the inter- 
mediated conclusion that: 

(44) Ind{7r+[/r7r+} = - J tro{(C/^ - J)[yr, Ur]}. 

The integrand of the last integral is finite, fact that can be seen from Propo- 
sition 3. Before we go further we will do two things. First, we will place an 
approximate projector squared vr^ in front of (C/p — /), which is legitimate 
since tTs acts as identity on the space where (U^ — I) is non-zero. It is easy to 
see that the conditions of the Proposition 4 are satisfied and for that reason 
we can move tt^ all the way to the right, inside tro. 

The second thing we do is, by using the Non-Commutative Residue Theo- 
rem, replace ?7p — / by ^{f{Ur) — f{I)), with f{z) analytic in a strip around 
the unit circle. This step will become crucial at a later point in the proof, 
but it must be done at this step. Thus, we arrived at 

(45) Ind{7r+t/r7r+} = -^/dF tro{(/(t/r) - /(/))[yr, C/rlvr^}. 
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U-32U3 .■■■--■■-■■■-■■■-■■■•■■■■^■■■-■■■■•^ 

u-22u2 - * * • * ^ 

u lu • 



E 

• ' • uiSu-i 

uaSu-B 



Figure 6. A graphical representation of ^inSt_„rt*n = 
^^(2n+ l)7rr(n). The top hnes represent the spectral repre- 
sentations of Ux„T,t_^ru% , which are discrete Heaviside func- 
tions shifted by n sites. All these operators act on 7^^• The 
sum of the top lines results in the stair like function shown 
at the bottom. 



Like in Ref. [5], we evaluate the commutator using Duhamel's identity 
[yr, Ur{Hr)] = - / dt^{t){l + it) dqx 

with (f) being the Laplace transform of e~'^^^F{x) _ which is a smooth 
function with support in [£'„ -|- 5,E^ — S]. Then: 

Ind{7r+f/rvr+} = dT f dt<^(t)(l + it) /J dqx 

tro{(/(t/r) - /(/))t/re-(i-'?)(i+^*)^r[yr,^r]7r,e-''(i+'*)^r7r,}. 

We will use now Proposition 4 to move e~'''^^^'^^^^^TTs all the way to the left, 
inside trg. Before doing that, we must be sure that the conditions of the 
Proposition 4 are satisfied. And indeed, the operators e~*(-'^+**)^i'7r57rr(0) 
and 7rr(0)e~''(^"'"**)^'"7rs are Hilbert-Schmidt, which follows from point (i) of 
Proposition 3 and the fact that tTs has same properties as Ur — I- For the 
same reasons, 

(48) 7rr(0)(/(t/r) - /(/))t/re-(i-'')(^+**)^^[yr, i^r]7r, 
and 

(49) ifiUr) ~ /(/))C/re-(i-'')(i+**)^^[yr,i^r]vr.^r(0) 

are Hilbert-Schmidt. After moving e^'^^^^^^^^^TTg all the way to the left, 
inside tro, we use the observation that all the operators to the left of [yr, Hr] 
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commute and that tTs can be replaced with the identity, to obtain 

Ind{7r+f/r7rit} = ^ J dT J dt^{t){l + it) dqx 
tro{(/(C/r) - /(I))^re-(i+'*)^^[yr,i^r]vr,} 

(50) 

= ^Jdr tro{(/([/r) - f{m'{Hr)[yT,HT]'Ks} 
= /dr tro{(/(i7r) - f{I))UrF'{Hr)[yT,Hr]'Ks} 

We now take 

(51) /(,) = _i_l_l,,>0, 
for which 6i=l. Then 

(52) Uf{Uv) - f{I))Ur = {Ur - I){Ur - (1 - e)J)-\ 
and by taking the Hmit e — ^ 0, we obtain: 

(53) Ind{7r+C/r7r+} = 27ri J dV tio{p{Hr)[yr , i?r]7r,}. 

The last equation is precisely the one of the theorem because tt^ acts as the 
identity over the space where p{Hr) is non-zero. 



5 Conclusions 



We presented a fairly elementary proof that the edge channels of Chern in- 
sulators are not destroyed by a rough edge. In the process, we demonstrated 
the use of the new topological invariant introduced in Ref. [3] for solving 
interesting problems in condensed matter. 

The technical estimates derived in this paper hold for very general tight- 
binding Hamiltonians, in particular for lattice systems with spin. Thus the 
paper provides the technical ground for applications to spin-Hall effect. In 
other words, the formalism can be applied line by line to this case too, 
the only thing that remains to be identified is what observable should one 
use instead of yr to get a non-zero value for the Index. This has been 
accomplished in Ref. [14]. 
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6 Appendix 

6.1 Proof of Proposition 1. 

This will be an elementary application of the Combes-Thomas principle [15]. 
We do not seek here optimal estimates, which could be obtained with the 
technique developed in Ref. [16], but rather the simplest proof of the above 
estimate. For this, consider the invertible transformation: 

(54) |na) ^ f/q|n,a) = e-'i"|n,a), q G R^. 
With the notation Hq = U-qHoUq, we have: 

(55) Hq= J2 [C^e'»™|n,a)(n + m,/3|+r™;e-'i'"|n + m,/?)(n,«|]. 

n,m,o,/3 

We write Hq = Hq + Wq, with 

^q= E [r™^[e'i--l]|n,a)(n + m,/3| 

(56) n,m,a,/3 

+r™;[e-'i--l]|n + m,/?)(n,«|]. 

Suppose we can show that Wq is small when q is small. To be more precise, 

assume that the operator norm ||W^q|| goes to zero as q(= jqj) goes to zero 
(recall that \\A\\ = sup where supremum goes over all ||/|| = = 

1). Then 

\\Hq-z\\ = \\Ho-z + Wq\\ 

(57) > \\Ho - z\\ - WWJ 

> dist{z,a{Ho)) - ||Wq||. 

The above inequality has substance only if the right hand side is positive, 
which is the case if q is small enough. This estimate gives 

"<^^-'"^ d.st(...W)-HI^..I| - 

which becomes useful in the following way: 

\{n,a\{Ho-z)-^\n',P)\e^(^-^') 
.^g. ={n,a\Uq{Ho-z)-'U.q\n',^)\ 

= \{n,a\{Hq-z)-^\n',P)\ 
<{distiz,aiHo))-\\Wq\\)-\ 
If we orient q parallel to n — n' we obtain 

(60) \{n,a\{Ho-zr'\n',(3)\< 



dist(z,cr(ifo)) - SUP|q|=g \\Wci\ 
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It remains to estimate ||VFq||. One can directly compute: 
(61) {n + m,P\W^\n,a) = (e"^- - l)[r^™ + r™^*]. 

We consider now two arbitrary unit vectors 

/ = En,a "^n," I «) i 9 = J2n,a ^n,a | n, «) , 

and we use the Schwarz inequality to derive an upper bound for ||M^q||: 

1/2 



(62) 



\{9mf)\ 



E6;.^(n',/3|T^q|/) 

n',/3 



< 



E|(n',/3|Wq|/)P 

n',/3 



(63) 



E 

n',/3 



X;(n',/3|Wq|n,a)an,a 



n,a 



1/2 



< 



1/2 



if(2L + l)2 X: EI(n',/3|Wq|n,a)|2|a„,aP ■ 

n',/3 n,a 

In the last step we used the fact that the number of non-zero terms in the 
sum over n and a cannot exceed K{2L+ 1)"^. We continue: 

1/2 

|2 



\{9mf)\ < 

(64) =(2L+l)v^ 



n,a m,/3 



EE(e-^'"-i)'|r^r + r2^;Pl«n,aP 



n,a m,/3 



1/2 



< (2L + l)x/Z sup 

a 

We can then take 



m* 1 2 



E(e-'i™-i)2|r^™ + r-^ 

m,/3 



1/2 



(65) C{q) = C^L + 1)Vk sup 

a,|q|=q 



5^(e-^™-i)2|r^™ + r™;p 



m,/3 



1/2 



which evidently decays to zero as g — > 0. 



6.2 Proof of Proposition 2. 



(i) First, we point out that 

(66) |(n,a|<^(i^r)|n',/?) = \{n, a\(P{H)\n' , (3) , 
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if both n and n' are in the "+" zone of the lattice. Thus we can work 
with H instead of Hr- Let us use the notation Rq{z)={Hq — z)~^ and 
R{z)={H — z)~^ . We make use of tlie following simple identity: 

(67) R{z) = Ro{z) + Roiz)AVRo{z) + Roiz)AVR{z)AVRo{z), 

and we do the functional calculus via the Stone's formula: 

[R{\ + ie) - R{X - ie)](j){X) — , 

R 

where the limit is in the weak operator topology. Using Eq. 67 and the fact 
that the support of ^ is in the gap of i^o, we obtain 

(n,a|</.(il)|n',/?) = lim E ! ^^ 

(69) (n,a|i2o(A)|m,a')(k,/?'|i?o(A)|n',/3)x 
(m, a'\AV[R{\ + ie) - R{\ - ie)]AV\k, /?'). 

It is important to notice that, due to the localization properties of AV, the 
sum over m and k can be restricted to the sites with first coordinate within 
the interval [-D-L+1,D+L-1]. We consider the following vectors: 

(70) Vm,a' = AF|m,a'), ^M' = Ay|k, /?'). 
The norm of the two vectors is bounded by 

(71) ||^™,„nMIV'M'|l <2if(2L + l)2max{|r™^|}. 

We will call Q the constant appearing to the right. With the notation 

(72) *m,a' = V'm,a'/IIV'm,a'||! ^k,/3' = V'k,/?'/ II V'k,/?' |l ) 

we continue as follows: 

{n,a\^{H)\n',(3) 



lim E IIV'm,a'|lllV'k,/3'||e-(^+«)(l™l+l"'-''l)/ 



(73) ^-*0m,k,a',/3' R 

.^(A)(n, a|i?o(A)|m, a')eMI"-"^l (k, f3'\Ro{\)\n' , /3)e(5+OI"'-k| x 

(*m,a'|[^(A + ie) - R{X - ^e)]|*k,/3'), 

where q>0 and ^>0 are chosen so that Proposition 1 applies with q replaced 
by q+^. We consider the following function, 

(74) K^^CW = <^(A)(n,a|i?o(A)|m,«')eMI"-H 

x(k,/3'|i?o(A)|n',/3)eMI"'-'*l 



2m 
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and observe that Eq. 73 can be written as: 

{n,a\(t>{H)\n',P)= ^ II Wll IIV'M'II 

m,k,a',/3' 

^^^^ xe-(i+Oi\^-^\+\n'-^)jF:^ll^^{X)di,{X) 

R 

where d//(A) is the spectral measure of H, projected on ^m,a' a^id ^k,/3') 
which are unit vectors. Thus, we can obtain an upper bound in the following 
way: 

(76) \{n,am)\n',P)\ < ^ e-^d"— l+l-'-^l) ^ sup\F^^SiX)\. 



m,k a',P' ^ 



Prom Proposition 1 we have: 



„ |paa'/3'/3/xN| . SUp^0(A) 

(77) sup|F_,„,(A)|<j^_^^^^^^j, 
Also, for the remaining sums we have 

(78) ^e-(«+^)l"-'"l < 5^e-«"i, ^ e-(«+^)l"'-^l < 5^e-«"i 

m k 

where is a positive parameter depending only on ^. If we put everything 
together we obtain: 

(79) \{n,am)\n',P)\ < Q^K'Slj^^^^^e'^^^^^^l 



(ii) Here we do the functional calculus using a technique introduced by 
Helffer and Sjostrand [17]. We consider the domain V in the complex plane 
defined by all those z with |Imz|<fo and Rez G [E^ + 5, E^ + 6]. We recall 
that the distance from the support of (j) and the spectrum of Hq is at least 
S. We will take vq less than this 6. Since (f)(e) is smooth, for any positive 
integer N, one can construct a function f:T>^C such that 

a) fN{z, z)=(l){z) when z is on the real axis. 

b) \dzfN{z,z) < aN\lmz\^. 

Such function is called an almost analytic extension of (p. Using such func- 
tion, one has 

(80) ^{H) = ^J d,fN{z,z){H - z)-'(fz. 

V 
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It is easy to see that the result stated in Proposition 1 applies equally well 
to the resolvent of H: 

(81) <".°i(^-^)->''^)^ dist(.r.w)''-cw ' 

for q small enough such that the denominator is positive. The difference 
between the two cases is that H may not have a spectral gap like Hq (see 
Fig. 5). Still, the result is of interest to us since it gives the behavior of R{X) 
for z in T> and away from the real axis. One can repeat the estimates on 
Wq, with the links crossing the contour T erased, and convince himself that 
the function ((q) remains the same. If z=u + iv, then dist{z, a{H))>\v\. 
Looking at the expression for (" function given in Eq. 65, we see that it 
behaves linearly with q for small values of q. Thus, if we take q=\v\/M, 
with M large enough, there is a > such that 

(82) dist(z, a{H)) - C(hl/M) > e\v\ > 0, 

for any z eV. This gives us an upper bound on the resolvent for all z in V: 

-|^;||n-n'|/M 

(83) {n,a\{H -u-iv)-^\n'p) < — . 

(/\v\ 

Then we can continue 

\{n,a\cj>{H)\n',(3)\ 

^Jdu J dv dzfN{z,z){n,a\{H -u-iv)-^\n',P) 

—vo 

vo 

<<^JduJdv vN-l^~v\n-n'\/M 

(84) fo 

^ ^ < a|e^o/M Jdu jdv ^iV-lg-^(|n-n'|+l)/M 

^ 

A,.iV , «o(|n-n'|+l)/M 

< gjyAM^ vo/M r J N-1 -X 







- 7r6»(|n-n'|+l)Jv'^ J ux x a - ^5)(|„_„/|+i)iv e 



(iii) Replace q by 2q at point (i) and and replace by 2A'' at point (ii) and 
take the product of Eqs. 11 and 12 to obtain: 

2 . ^ „ g-2g(ni+n'i) 



(85) |(n, a\ct^{H+)\n', P)f < A{2q)B2N _ ^ ^^2iv • 

The statement follows because we can replace |n — n'| by |n2 — any time. 
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6.3 Proof of Proposition 3. 



Proof, (i) We will show that the Hilbert-Schmidt norms are uniformly 

bounded: 

(86) \\{Ur - I)Mn)\\is = TV{7rr(n)([/^ - I){Ur - I)Mn)}. 
Since 

(87) {U^ - I){Ur -I) = {e'^'''^^^^^ - i)(e-^^'^^^^^ - 1) 

and the function (e27rjF(a;) _ -|^-j^g-27riF(a;) _ 2^ smooth and with support in 
the interval [E^ + S, E^ — S], we can apply Proposition 2 to conclude at this 
step that 

(88) {n,a\{U^ - I){Ur - I)\n',(3) < ^(g)e-«("i+"i). 
Then 

TV{7rr(n)(C/i^ - I){Ur - I)Mn)} 
= E {n,a\{U^-I){Ur-I)\n,a) 

ni>j„,a 

(89) <K Yl ^(g)e-2«"i 

ni>7„ 

<K Yl ^(g)e-2«"i < oo. 

ni>~D 



(ii) 
(90) 



\\Mn)iUr - I)Mn')\\is 
= Tr{7rr(n')(t/? - I)Mn){Ur - I)Mn')} 
E E E\{{ni,n),a\{Ur-I)\{n[,n'),P)\'' 



ni>jn ni>7„/ a,l3 

and by applying Proposition 2 point (iii) we have 

||7rr(n)(C/r-/)7rr(n')|||s 



(91) ni>7nn'i>7, 



<P=,fti)w(^^E_^^iv(5)e-^-j 



2 



(iii) We proceed as follows. 

IIPr,^r]||L= E lkr(n)[Sr,C/rKr(nO|||s 

(92) 



n.n' 



E \\Mn)iUu,-I)Mn')fus= E Krin,n'). 

n-n'<0 n-n'<0 
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In the last two sums, we must exclude the term n=n'=0. If wc take N > 3 
at point (ii), the final sum is convergent and uniformly bounded. 

(iv) We use the following equivalent expression for the commutator [yr, Ur]- 
(93) [yr, Ur] = J2in- n')7rr(n)(C/r - I)T^r{n') 

n,n' 

to obtain: 

\\[yT,Ur]Mn)\\is = T.\n-nnMnm-I)Mn')\\is 

n' 

= ^ |n — n'\^KY{n, n'). 

n' 

If we take A?^ > 3 at point (ii), the final sum is convergent and uniformly 
bounded. 

6.4 Proof of Proposition 4. 

M 

Proof. Let 7r(M) = X) 7i"r(n). Then 7rr(0)^r7r(M)Br7rr(0) are trace class 

-M 

and 

/ dV Tr{7rr(0)Ar7r(M)5r7rr(0)} 

M 

= JdT E Tr{7rr(0)^r7rr(n)i?rvrr(0)} 

(y4j n=-M 

M 

= Jdr Y: 'IV{7rr(n)Sr7rr(0)^rvrr(n)} 

n=-M 

At this point we use the invariance of the trace (on trace class operators) 
under the unitary transformations to continue: 

M 

... = JdF J2 Tr{unMn)BTMO)ArMn)<} 

n=-M 

M 

(95) =IdT Y TV{7rt„r(0)Si„r7rt„r(n)^t„r7rt„r(0)} 

n=-M 

M 

= Jdr E TV{7rr(0)5r7rr(n)^r7rr(0)}. 

n=-M 

At the end of above argument we used the invariance of dT relative to the 
transformations tn- We then have that 

JdT Tr{7rr(0)^rvr(M)Br^r(0)} 

(96) 

= JdT TV{7rr(0)Sr7r(M)^rvrr(0)}, 
and the affirmation follows by letting M go to infinity. 
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